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1. Introduction 



One of the outstanding open challenges for theoretical fluid mechanics in the 21^* century is to derive rigorous 
results for turbulence directly from the fundamental equations of motion, the Navier-Stokes equations, 
without imposing ad hoc assumptions or uncontrolled closures. Exact results are extremely rare, but it 
is possible to derive rigorous and physically meaningful limits on some of the fundamental physical variables 
quantifying turbulent dynamics and transport. The bulk rate of energy dissipation is one such quantity of 
particular interest due to its production as a result of the turbulent cascade in the high Reynolds number 
vanishing viscosity limit. The derivation of mathematically rigorous bounds on the energy dissipation rate, 
and hence also a variety important quantities such as turbulent drag coefficients and heat and mass transport 
rates, has been a lively area of research in recent decades. 

Beginning in the early 1960s, L.N. Howard and F.H. Busse pioneered the application of variational 
approaches for the derivation of rigorous — and physically relevant — bounds on the dissipation rate for 
boundary-driven flows; see their reviews [201 In the 1990s, P. Constantin and the senior author of 
this paper introduced the the so-called background flow method [HI E] based on an old idea by Hopf 19 . 
The background method was soon improved by Nicodemus et al '22' who introduced an additional variational 
'balance' parameter, and by the late 1990s Kerswell j21j had shown that the background method equipped 
with the balance parameter is dual to the Howard-Busse variational approach. Those theoretical techniques 
have been applied to many flows driven by boundary conditions, including shear flows and a variety of 
thermal convection problems pi H [111151 [TTl 1^ 1^ . 

Attention has recently turned as well to the derivation of quantitative variational bounds on the energy 
dissipation rate for body-forced flows. In these systems, the bulk (space and time averaged) dissipation 
rate per unit mass e is proportional to the power required to maintain a statistically steady turbulent state. 
While body forces may be difflcult to realize in experiments, they are easily implemented computationally 
and are the standard method of driving for direct numerical simulations (DNS) of turbulent flows. 

Childress et al 2 applied a background-type method to body-forced flows in a periodic domain, focusing 
on dissipation estimates in terms of the magnitude of the applied force. In dimensionless variables they 
bounded e in units of {F^£y/^, where F is the amplitude of the applied force per unit mass and i is the 
(lowest) length scale in the force. The estimates were given in terms of the natural dimensionless control 
parameter, the Grashof number, Gr :— FH? jv^^ where v is the kinematic viscosity. In practice, e is often 
measured in inviscid units of U'^ j i as a function of the Reynolds number i?e = UH/v, where J7 is a relevant 
velocity scale — an emergent quantity when the force is specified a priori. In both cases the dissipation is 
bounded on one side by that of the associated Stokes flow 17 . When bounds are expressed in terms of Gr, 
the Stokes limit is an upper bound, whereas when the estimates are in terms of Re it is the lower limit. 

Foias ^31 was the first to derive an upper bound on 

in terms i?e, but with an inappropriate prefactor dependence on the aspect ratio a = where L is the 
system volume, generally an independent variable from £ (see also jl5l 116 ) ). That analysis was recently 
refined by Foias and one of the authors of this paper |^ to an upper estimate of the form 

/? < ci + ^ , 
He 

where the coefficients Ci and C2 are independent of F, £, v and a, depending only on the "shape" of the 
(square integrable) body force. (This in consistent with much of the conventional wisdom about the cascade 
in homogeneous isotropic turbulence theory ^1 IIUL ITl] as well as with wind tunnel measurements ^] 
and DNS data (28j.) Most recently, that approach was developed further by deriving a mini-max variational 
problem on the time averaged dissipation rate for a particular domain geometry [S]. Moreover, the variational 
problem was solved exactly at high Reynolds numbers to produce estimates on the asymptotic behavior of 
the energy dissipation as a function of Re including the optimal prefactor. 

In this paper we extend the results in |H] by introducing a balance parameter c, the analog of the 
variational parameter introduced by Nicodemus et al 1 22112311^ for the background method. This parameter 
controls a balance between the quantity being bounded, the manifestly positive definite energy dissipation 
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rate proportional to the norm of the rate of strain tensor, and the indefinite quantity derived from 
the power balance that is ultimately being extremized. Specifically we consider the flow of a viscous 
incompressible fluid bounded by two parallel planes with free-slip boundary conditions at the walls and 
periodic boundary conditions in the other two directions. The flow is maintained by a time-independent 
body force in the direction parallel to the walls. First we derive the Euler-Lagrange equations in the case 
c — Q (where the variational principle coincides with the one in 0) and solve them numerically at finite Re. 
The full (c > 0) Euler-Lagrange equations are quite complicated but they can also be solved numerically by 
using Newton method with the c = solution as an initial guess. 

The rest of this paper is organized as follows. In Section[21we introduce the problem and its variational 
formulation following In Section |3| we present the augmented variational problem and derive the 

variational equations, explaining how we go about solving them. In Section ^ we collect our numerical 
results, and in Section we summarize the results discussing the challenges of this approach and future 
directions for research. 



2. Statement of the problem 

2.1. Notation 

Consider a viscous incompressible Newtonian fluid moving between two parallel planes located at y = and 
y = i. Denote x the stream-wise direction and z be the span- wise direction. The velocity vector field satisfies 
free-slip boundary conditions at the two planes bounding the flow. We impose periodic boundary conditions 
in the other two directions. The motion of the fluid is induced by a steady body force f along the x axis 
varying only in the y direction. 

The motion of the fluid is governed by Navier-Stokes equation 

^ + (u.V)u + Vp^i^Au + f (1) 

and the incompressibility condition, 

V • u = . (2) 

Here p(x, t) is the pressure field, and Re := '^'""^ is the Reynolds number, where C/rms is the root-mean 
square velocity of the fluid. The problem is non-dimcnsionalized by choosing the unit of length to be i and 
the unit for time to be £/J7ims- Let (•) stand for the space-time average. With this choice of units the velocity 
of the fluid u(x, t) = {u, V, w) is space-time L^-normalized to 1: 

{\u\^) = {u^ + + w^) = I . (3) 

Given e, is the space-time average energy dissipation rate in physical units, the non-dimensional energy 
dissipation rate P is defined 

(4) 

yy rnis 

The body force f in (^) has the form 
f(x) = F0(2;)e, , 

where the dimensionless shape function ^ : [0, f ] — > M has zero mean and satisfies homogeneous Neumann 

boundary conditions, and is L^-normalized: 

.1 .1 

0(y)dy = O, 0'(O) = O = 0'(1) , / c^{yfdy = l. 

Jo 

Now let $ e H^{[0,1]) (where iJ''([0, 1]) is the space of functions defined on [0,1] with L^-integrable p*'* 
derivatives) be the potential defined by 

$' = -0 , $(0) = $(1) . 

(Note that we are free to impose homogeneous Dirichlet conditions on <& at both boundaries due to the zero 
mean condition on <p.) 
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The spatial domain is {x, y, z) E [0, L^] x [0, 1] x [0, L^] where and are the (non-dimensionaHzed) 
lengths in x and z directions. Free-slip boundary conditions at the walls are realized by 

. = 0, 1 = 0^1^ at, .0,1. (5) 

2.2. Variational problem for the energy dissipation rate 

Here we follow jBj to derive the variational problem for upper bounds on the energy dissipation. Multiplying 
Navier-Stokes equation 0J by u, integrate over the spatial domain, and average over time to obtain the 
energy dissipation rate 

/? := -^(|Vu|2) = (f • u) = F{c^u) = -F{<i>'u) . (6) 
Re 

To remove the explicit appearance of the amplitude F of the body force, multiply by a vector field 
of the form ip{y)ex, where the multiplier function tp € H'^{[0, 1]) satisfies homogeneous Neumann boundary 
conditions V''(0) = = and is not orthogonal to the shape function cf). That is, {(pip) 0. We will also 

use the derivative of -0 

^ = ip' e H\[0,l]) 

which satisfies homogeneous Dirichlet boundary conditions ^(0) = = ^'(1) and is not orthogonal to the 
shape potential i.e., ($^) — ^ 0. We will call ^' a test function. Take the scalar product of ^ with 
ip{y)ex, integrate over the volume (integrating by parts by utilizing the boundary conditions) and take the 
long-time average to see that 

- = —(^'u) + . (7) 

Re 

Express the amplitude F of the body force from Q) and insert into the expression for the energy dissipation 
© to obtain 

^= W) ■ 

2.3. Mini-max upper bounds for [3 

A variational bound on (3 may be obtained by first maximizing the right-hand side of ^ over all unit- 
normalized divergence- free vector fields u that satisfy the boundary conditions |SJl, and then minimizing 
over all choices of test functions 4* € H^{[Q, 1]) satisfying homogeneous Dirichlet boundary conditions. Then 
any solution of Navier-Stokes equation will have energy dissipation rate /3 bounded from above by 

Oftlite) = mmmax ■ . (9) 

* u ($11/) ^ ' 

In order to study the bound above, the authors of 8 first evaluated (exactly) 

. {"^uv) 
Pb[oo) :=mmmax — — , 

* u (<I)\p) 

and then used this result to analyze the behavior of (3b{Re) for finite Re. Since we are going to generalize 
that approach, we briefly recall the analysis: 
The evaluation began with the proof that 

max (<i>'u) (*to) = sup |*(y)| • (10) 

V27 ye[o,i] 

This was accomplished by showing that the right-hand side of (|10|) is an upper bound for ($'w)(^'uu) for 
any u in the class of vector field considered, and then explicitly constructing a sequence of unit-normalized 
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divergence- free vector fields u*^*^^ = {u^''\v^''\w^^'>) satisfying the boundary conditions ^ such that u*^*^) 
saturate this bound in the limit k — > oo, i.e., 

lim (^'m^*^)) = sup \^{y)\ . 

fc^oo v27 j/e[o,i] 

The precise form of u'^'^'^ is 

u'^''Hy,z) = .9fe(y)\/2sinfcz- -i=$'(y) 

v'^''Hy,z) =gk{y)V2 sin kz (11) 

w^''\y,z) = j:9'k{y)V2 cos kz , 

where the sequence consists of smooth functions approximating as /c ^ oo a Dirac 5 function with support 
centered at the points where the function ^' G H^{[0, 1]) reaches an extremum, and normalized as 

2 1 /2\ 1 

Note that the function ^' e -^^([0, 1]) is continuous and hence it reaches its extremum in [0, 1]. Moreover, 
since ^'(0) = = ^(1) and at the same time 5* is not identically zero, a point where reaches an extremum 
must be in the open interval (0, 1). 

Following (|10|) . it was proved that if $ e ^^([0, 1]) changes sign only finitely many times, then 

Pbioo) = —= min sup ' - 



which is achieved for the choice of test function ^ = sign<I>. While sign<i> is not in if^([0, 1]), it can 
be approximated arbitrarily closely (in the sense of pointwise convergence) by a sequence of functions in 
i?i([0,l]). 

In [21, the authors considered test functions which are "linearly mollified" approximations of sign$, 
i.e., continuous piecewise linear functions approximating sign<i> by replacing the jumps of sign$ by lines of 
slope ±^ connecting the values —1 and 1 (see Figure 1 in |H]). Finally, for finite i?e, it was shown in |Hj that 
by choosing 8 ^ 0(Re~^^^), the dissipation rate for $ £ -^"'^([0, 1]) behaves for large Re as 

Pb{Re) <l3b{oo)+OiRe-^/^) . 

If $ is smooth (i.e., <I> has a bounded derivative and so behaves linearly around its zeroes), then by taking 
i5 ~ 0{Re^^^^) it was shown as well that 

l3b{Re)<l3b{^) + 0{Re-^'^) . 

3. Improved variational principle 

3.1. Introducing the balance parameter 

Let c e [0, oo) be arbitrary. Multiply © by 1 + c and add it to /9 = -^(iVup) multiplied by — c. The result 
is 

/?=(!+ -y-^^^^f^ ]^(ivun . (12) 

Now we will obtain bounds on the energy dissipation by applying a mini-max procedure to the functional in 
the right-hand side above. 

The parameter c provides more constraint on the variational procedure than the case considered in 
The space-time average of |Vup is multiplied by — c < so that for a velocity field with a large gradient 
(like the one of the form Hll|l when g^ tends to a Dirac 5 function), the right-hand side of (|12|l will become 
smaller. 
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While performing the maximization procedure we have to incorporate two exphcit constraints on the 
velocity vector fields: the unit-norm condition ((JJ and incompressibility The former one is easy to 
implement by adding a term with Lagrange multiplier A which is a number (i.e., does not depend on x and 
t). Incompressibility, however, requires introducing a Lagrange multiplier (a "pressure") that is a pointwise 
function which makes the variational problem very difficult to analyze. So instead we will restrict the class 
of velocity fields u over which we maximize to fields that are automatically divergence-free. 

The functional incorporating the normalization constraint is 



L[u] := (1 + c) 



-1) 



(13) 



The class of velocity fields u we will consider is a generalization of (jllfl : 
u(jj, z) — U (y) V2 sin kz + A(y) 

v{y,z) ^V{y)V2sinkz (14) 

w(y, z) = Y y'{y) V2 cos kz , 
k 

where the functions U, V , and A satisfy the boundary conditions 

C/'(a) = F(a) = F"(a) = A'(a) = , a = 0, 1 . (15) 
Note that the vector field u defined in (|14|) is automatically divergence- free. 

This class of velocity fields u (|14|) is restrictive, but in our opinion it constitutes a physically reasonable 
ansatz. It has been observed for plane parallel shear flows that the first modes to lose absolute stability have 
only cross-stream and span-wise variation with no dependence on the stream-wise coordinate x. Moreover, 
the parameter k in p4(l can take any real value, so this does not impose any restriction on the wavelength of 
the pattern in span-wise (z) direction. Note also that the case of very high Reynolds numbers corresponds 
to the choice c = (see (|13|l '). and in this case the family H14|l will tend to the family 111|) which we know 
achieves the upper bound on the dissipation at infinite Re. All these considerations make the choice of 
the family (|14ll quite reasonable. In the spirit of full disclosure, however, we reiterate emphatically the 
assumption that we make in the analysis that follows: 

Ansatz: We assume that the maximizing vector fields for the functional (|13|l have the functional form (|14|l . 

In terms of U, V, and A, the expression H12I) for the energy dissipation reads 

($'A)(*C/y + -^*'A) 



/3[C/,y,A] = (1 + c)- 



- (k'^U^ + k^V^ + U'^ + 2V'^ + 
Re \ k^ 

and the functional L[u] 113|l taking into account the normalization constraint becomes 

A 

2 

The Euler-Lagrange equations for U, V, A are 



L[U, V, A] = f3[U, V,A] + ^{U^ + V^ 



Re 



= 



(16a) 



2c 
' Rek^ 



V" 



4c 
Re 



A 

fc2 



V" 



A-^V+(i+c)is:#w 



Re 



= 



(16&) 



2c 
Re 



1 



A" + AA + — (1 + c) 
Re 



- c]-, f - 



AWV) 1 ,(*'A> 



$' = 



where the "eigenvalue" A is to be adjusted so that the triple {U, V, A) satisfies the normalization 
+ + -^V'^ + A^) = 1 



(16c) 



(17) 
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3.2. Exact solution at finite Re for the case c — 

In the case c — 0, the Euler-Lagrange equations (|16a|) . I|166|) . HKicH become 
($'A> 



XU- 

_x_ 

AA- 



V" + XV 







($'A) 



1 (*'A),x. 



Re 



(^'C/y> 1 (*'A) 



$' = 



(18a) 
(18&) 
(18c) 



Then the equations for U and A are algebraic equations, so the only boundary conditions that have to be 
satisfied are 

V{a) = for a = 0, 1 . (19) 

We can solve the boundary value problem (|18af) . Ijl86|) . (|18c|) . (|19() explicitly. First, expressing U from 
()18a|) . and substituting into Ijl86|) . we obtain the following boundary value problem for V: 

^2 



fc2 



V = E^ 



■V , V{0) = V{1) = 0, 



where we have set 
E 



($'A)y(¥^ 



(20) 



(21) 



For each choice of test function 5" we obtain a sequence of functions Vn and numbers En, n — 1,2,3,.... 
For each n, the numbers £"„ and the functions Vn depend on Re, k, and the choice of test function 5*. The 
functions A„ are (see the Appendix for a derivation) 



An(y) 



1 

71 



'1 



£■2 



and the functions C/„ are 

Un{y) = - 



■■Vn{y) 



3i?ey(*2) 



En 



Re^/Jp 



:*'(y) , 



(22) 



(23) 



In the derivation of (|22|l we used the normalization condition (|17|l so that it is automatically satisfied. Then 
the (non-dimensional) energy dissipation rate is 



Pn = 



3($*) ^ ^ 3($*)(*2) 
1 



1 



1 



^-2 



/2\ 
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($'r)2 



3/2 



(24) 



What remains to be done for a given shape potential and multiplier function is to find the solutions for V 
and E. This we do numerically. 



3.3. Finding the velocity profile and energy dissipation for c > 

Suppose that we have found the functions Un^ (|23|l . (I2UII . and Kn^ (|22|) satisfying the Euler-Lagrange 
equations (|18af) . (jl8fe|) . (|18cf) and the boundary conditions 119|l in the case c = 0. In order to find the solution 
Un, Vn, An of the boundary value problem 116all . 1166II . 116cll . H15|l that satisfy the normalization condition 
(I17II for c > 0, we use Newton method with C/,l"\ vJi"\ A'^^^ as initial guess. 

According to the general methodology of the mini-max procedure, we have to first maximize the 
expression for the energy dissipation rate /3 over all allowed velocity fields u ()14|l. and then to minimize 
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maxu (3 over all allowed functions With our ansatz for the form of u, maximizing over u means maximizing 
over all real values of k. Then having found the maximum of /? over fc, we minimize over both ^' and the 
balance parameter c > 0. In practice we have to choose a particular family of test functions ^' depending on 
a small number of parameters, and minimize over those parameters and c. We will take a 1-parameter family 
of test functions 'i's (given explicitly in (|27|l below) where the parameter (5 is a measure of the thickness of 
a "boundary layer" . 

Let /3(i?e, i5, c, k) be the mini-max upper bound for the turbulent energy dissipation as a function of the 
Reynolds number Re, the parameter S of the family '^g, the balance parameter c, and the wavenumber k. 
Define (3* {Re, 5, c) to be the maximum over k of (3{Re, S, c, k), and k*{Re, 5, c) to be the value of k for which 
(3{Re,5,c,k) attains this maximum. Then 

P*{Re,S,c) :— ma.x (3 {Re, S,c,k) , k*{Re,S,c) := argmax /3(i?e, <5, c, •) . (25) 

k 

After maximizing over k, i.e., over the family of velocity fields u (I14|l . we minimize over the parameter S of 
the family of test functions '^s, and the balance parameter c. That is, we compute 

^b{Re) -.^ mm P* {Re, S,c) , {S* {Re), c* {Re)) := urgmm f3* {Re, ■,■) . (26) 

(5,c) 

4. Numerical results 

4-1. Numerical example and implementation 

As a specific model to analyze we chose the same shape function as in (Hj: 
$(?/) ~ — sin Try , (/)(j/) — —^'{y) — —V^cosny . 

TT 

In the test functions 'fa were chosen piecewise linear but not continuously differentiable. For 
computational reasons we replace them with the smooth family 

^^>s{y)^{^-c-y/^){^-e-<^^-y'>^^) , s>o. (27) 

The functions H27|) satisfy the boundary conditions '^^(O) = = ^(s(l)- 

The boundary conditions of the Euler-Lagrange equations naturally suggest the use of Chebyshev 
polynomials as interpolants to implement a pseudo-spectral scheme [23 to solve these equations. The 
Matlab differentiation matrix suite '30' simplifies the implementation by providing routines to discretize and 
represent differentiation operators as matrices. Differentiation of a function then becomes multiplication of 
the differentiation matrix with the vector of the function values at those Chebyshev nodes. However, the 
discretized equations are still nonlinear in the c 7^ case. We started with the c = equations which are 
solvable as a linear eigenvalue problem H20|) . Then the standard Newton's method was applied to these 
solutions and iterated to solve the nonlinear equations pfiojl . Ijl6&|l . (|16c|l . The Jacobian matrices needed in 
the Newton's method were computed by a simple forward difference scheme. Throughout all computations, 
128 and 64 Chebyshev nodes were used (the differences between the results for these choices of number of 
nodes did not exceed 10^^). 

To illustrate the typical geometry of the flow, in Figures^andl^l we show the three coordinate projections 
and the 3-dimensional view of typical integral lines (i.e., solutions of {x,y,z) = {u,v,w) for {u,v,w) given 
by H14|) ') of the maximizing flow field for Re = 50 and Re = 1000, respectively. The values of the parameters 

5, c, k, for the fields shown are the ones that give the optimal bound, j3b{Re) given by H26|) . 

As an example of the mini-max procedure, we show in Figure 13 the upper bound on the dissipation for 
Re = 50 obtained by using as a test function from (|27|l with 5 = 0.04; the bound is given as a function 
ofce [0,1] and fee (0,25]. 

In Figure^ we show the bound on the dissipation /3 for Re = 50 as a function of the balance parameter 
c for different values of the span- wise wavenumber k; the data presented have been obtained with 'i's with 
S = 0.04. The figure illustrates the general behavior of /? as a function of k and c - namely, for small k, the 
value of P increases with c, while for larger k, f3 decreases with c. Clearly, the family of lines in the figure has 
an envelope - this envelope is the graph of the function /3*(50, 0.04, c) (|25|l . Having obtained the envelope, 
we find the minimum value of /3*(50, 0.04, c) - this is the mini-max value we are looking for; this point is 
labeled with Pi, in Figures 01 (where it is the saddle point) and^ 
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4-2. Results 

In Figure we present the bounds from previous papers, as well as our new numerical results. The dotted 
straight line represents the lower limit on the dissipation corresponding to Stokes (laminar) flow, 

PStokcs > -TT- ■ 

Re 

The dot-dashed line in the upper part of the figure is the bound following "5 for this problem obtained with 
^ = $: 

The thin solid line shows the "non-optimal" bound from |S] (equation (3.14) in |S]), 

PDES, non — optimal — 7== , . 5 

V27 Re 

while the long-dashed one gives their "optimal" estimate (obtained from equation (3.12) in [5] by first 
minimizing over ^ and then plugging ^' = $): 



V2Tr f X 3/2 ^2 / 

/5dES, optimal < —7= 1 + TT^Jl + 1 



(Note that this line bifurcates from the lower Stokes bound at Re = ^/2tt « 4.4429). The thick solid line 
starting from Re w 178 is the best upper bound for high values of Re from Theorem 1 of |8j: 

TT^ 5(67r2)i/5 3.73089 , , 

/?DES, Ti™. 1 < ^ + 44/5 R,./, - 0-67154 + . (28) 
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Figure 4. Bound on /3 for Re = 50 (obtained with ^0.04) as a function of c for several values of k. 



The circles in the figure give our new numerically determined upper bounds on f3 with the choice 5* 
the crosses represent our numerical results for the choice (|27|) of ^'5. 



$ and 



In Figure we have plotted /3 



(circles), k* (stars), and c* (x's), versus Re for the values of the 



dissipation bound obtained using the function '^s from H27|l . We see that k* 



'Re, 



Re 



and from the 



figure we observe that (3 



/216 



also behaves like a power of Re. In the figure we illustrate these behaviors 



by showing the straight lines 

-2 2.158 



k* = 1.0 /Re 



5.0 

'Re 



5. Concluding remarks 

We have derived new bounds on the energy dissipation rate for an example of body-force driven flow in 
a slippery channel. The fundamental improvement over previous results came from the application of the 
balance parameter in the variational formulation of the bounds, together with numerical solution of the 
Euler-Lagrange equations for the best estimate. 

In Figure the results of this analysis are compared with the direct numerical simulations of the 
three-dimensional Navier-Stokes equations first reported in [Sj- Over the Reynolds number range 100-1000 
where the data lie, the best bounds derived here, using the balance parameter and minimization over the 
(restricted) family of multiplier functions result in a quantitative improvement over the previous rigorous 
estimates. We observe that the measured dissipation is a factor of 3 to 4 below the bound, which should be 
considered nontrivial given the a priori nature of the estimates derived here. Presumably a full optimization 
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Figure 5. Upper and lower bounds on /3. 

over possible multiplier functions 5* would result in a further lowering of the estimate at lower values of 
Re, producing a bound that intersects the lower Stokes bound right at the energy stability limit (which we 
compute to be at Re — 2tt). We note from Figure [S] that the bounds computed with $5 tend to agree with 
those computed using ^ — "if a,t lower values of Re, indicating that both trial functions are about the same 
"distance" from the true optimal multiplier. 

At higher Reynolds numbers the optimal solutions computed here converge rapidly to the asymptotic 
bound Pb{oo) computed analytically in [5]. Indeed, the bound derived here approaches the asymptotic limit 
with a difference vanishing Re~^ '^^. This particular scaling of the approach to the asymptotic limit helps 
to understand the role that the balance parameter plays to lower the bound: while a naive estimate suggests 
that the approach might be C'(i?e~^), the faster convergence may be attributed to the interplay of the 
c ~ i?e~^ and k V Re scaling in the prefactor and the subtracted term in (|12|l . 

There are several directions in which this line of research could be continued. One is to develop more 
reliable and accurate analytical methods for estimating the best bounds at finite Re. This would probably 
involve asymptotic approximations for small but finite values of Re~^ which could lead to more general 
applications for other variational problems as well. Another direction would be to develop methods to 
determine the true optimal multiplier function at finite Re. The motivation there would largely be as a 
point of principle, to demonstrate that the full min-max procedure can indeed be carried out — at least for 
simple set-ups such as those considered here. Finally, going beyond the simple sin Try forcing considered in 
this paper there remains the question, first posed in [Sj, as to the connection between the optimal multiplier 
and the true mean profile realized in direct numerical simulations. Specifically, the question is whether there 
is a sensible correspondence between the shape of the optimal multiplier and the mean profile for general 
force shapes. The idea is that the optimal multiplier contains information about the extreme fluctuations 
that might be realized in a turbulent flow, and some of those features may correlate with the statistical 
properties of the flows. 
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Appendix: Derivation of the expression ()22|) for A 

In this Appendix we show how to derive the expression (|22|) for A in the case c = 0. First exclude U from 
(fTFcjl with the help of ((TFa|) : 



Now multiply the equation for U H18a|l by —U, add it to the equation for V (I18b|) multiplied by V, 
and integrate the resulting identity to get the equidistribution property {U^) = {V^ + -p-l^'^), so that the 
normalization condition (|17l) can now be written as 




(Al) 




(A.2) 



Multiplying ()20|l by V and integrating using the boundary conditions Ijiyf) . we obtain 
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Figure 7. Comparison between theoretical results and DNS data (same symbols as in Figure IBt. 



which, together with the new normahzation l|A.2p . yields 
2 " (*2) • 

This expression and the definition of E H21|) allow us to write the coefficient of the term of order Re'^ in the 
right-hand side of IjA.ip as 

($'A)(^V^) _ E^i^^V^) _ 1 - (A^) 
A2($5')2 " (5'2)($'A) ^ 2($'A) ■ 
Using the above relationship and expressing the Lagrange multiplier A from 12111 . we can rewrite (|A.1|1 as 

. /l-(A2) 1 E (*'A)\ , IE. 

A = ^ ; ^ '- ; . (A. 3) 

2($'A) Re ($'A) J Re ^ ^ 

Let ^ be the coefficient of $' in (|A.3|I . i.e., 

A:=A^<i>'-^-^*'. (A.4) 

Re 

From this expression we easily obtain (recall that ($'2) = 1) 



($'A) = - 
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Plugging these expressions in the definition of the coefficient /i, 
_ 1 - (A^) 1 E i^'A) 
^ ^ 2($'A> " ^ (^-'A) ' 

we obtain the following quadratic equation for 
The "physical" solution of this equation (the one that has the right behavior in the limit Re — s- oo) is 



Plugging this into IIA.4p . we obtain the desired expression H22|) . 
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